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Abstract 

We investigate the properties of the Sakai-Sugimoto model at finite magnetic field and baryon 
chemical potentials. We show that in a finite magnetic field, there exists a spatially homogeneous 
configuration carrying finite baryon number density. At low magnetic field and baryon chemical 
potential the equation of state of the matter coincides with that obtained from the chiral perturba- 
tion theory Lagrangian with an anomalous term. We discuss the behavior of the system at larger 
magnetic fields. 
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I. INTRODUCTION 

Recently, gauge/gravity duality has been used extensively to investigate properties 

of strongly coupled gauge theories at finite temperature and/or density. The popularity of 
this method is due to its ability to calculate in the strong coupling regime. The main 
disadvantage is that the models that can be solved, for example the Af = 4 super Yang- 
Mills (SYM) theory, do not coincide with QCD, typically containing additional degrees of 
freedoms like adjoint fermions or scalars. 

While finite-temperature A/" = 4 SYM plasma has many features reminiscent of QCD 
plasma at temperatures not too large compared to the deconfinement temperature, it is more 
difficult to construct a holographic model of cold nuclear or quark matter. One problem is 
that at large Nc nuclear matter is a crystal instead of a liquid. This fact finds refiection in 
the Sakai-Sugimoto model [J], where baryons are 5D instanton particles and nuclear matter 
is a crystal of such instantons [sl], which is necessarily inhomogeneous 6]. 

In this paper, we investigate the possible gravity dual of magnetized nuclear matter. In a 
recent study [tI , it was found that in the chiral limit of massless quarks, at any magnetic field 
the ground state of finite-density matter is not a crystal, but a spatially homogeneous phase. 
At low density such a phase is characterized by a finite gradient of the vr^ field, Vvr^ ^ [in 
the parametrization where the chiral condensate is proportional to exp(z7r''r'^//^)]. At finite 
quark masses this state becomes a stack of vr^ domain walls. The 7r° domain wall is locally 
stable in magnetic fields stronger than Bq, and is energetically more favorable than nuclear 
matter above a magnetic field Bi, where both Bq and Bi vanish in the chiral limit. 

The treatment of Ref. [3] relies on the use of chiral perturbation theory, including the 
appropriate WZW term in the presence of an electromagnetic field and a baryon number 
chemical potential. It is valid only for sufficiently small magnetic field and baryon chemical 
potential. 

In this paper, we search for a similar solution in the Sakai-Sugimoto model of holographic 
QCD [4]. The Sakai-Sugimoto model is an application of the AdS/CFT conjecture involving 
a system of Nf D8-D8 probe brane pairs in a D4 brane background of type IIA string 
theory. The model exhibits chiral symmetry breaking and thereby reproduces much of the 
low-energy physics of massless QCD, such as the octet of pseudoscalar Nambu-Goldstone 
bosons. Including the Chern-Simons (CS) term of the probe brane action is equivalent to 



2 



including the effects of the axial anomaly on the field theory side. While the model has 
been used to investigate properties of the vacuum and the thermal state (with zero chemical 

n Q 

potentials) in an external mag netic field |8|, M\, the case when both the magnetic field and 
the baryon chemical potentials are nonzero has not been considered. 

Because the Sakai-Sugimoto model incorporates the axial anomaly into a theory of mass- 
less pions, one would expect to find, at least at small magnetic fields and baryon chemical 
potentials, a solution similar to the one found in [3] . The purpose of this paper is to demon- 
strate that solution. 

We found, as expected, that at low B the results of Ref. [7| are reproduced. At larger 
values of the magnetic field, the quadratic approximation to the DBI action of the probe 
branes can no longer be trusted. However, we can still consider the quadratic action at large 



-B as a bottom-up AdS/QCD theory [lO|, lUl, ll2| . In this sense, we unexpectedly discovered 



that in the opposite limit of large B the zero-temperature thermodynamics of matter with 
finite baryon density is identical to the thermodyanamics of free quarks, which fill energy 
levels in the lowest Landau level. It is rather surprising given that the supergravity limit 
corresponds to the strong coupling regime in field theory. 

In Sec. [ni the necessary facets of the Sakai-Sugimoto model are reviewed. In Sec. Illlt the 
domain wall solution is presented. We give interpretations of the result in Sec. [IVl Section IVl 
provides a concluding discussion. 

Note added: While this paper was being completed, we learned that a similar investigation 
was also being carried out in [l^, which contains some overlapping results. 

II. REVIEW OF THE SAKAI-SUGIMOTO MODEL 



In iJ], a way to holographically model nonsupersymmetric pure Yang-Mills theory was 
presented. Although the initial motivation of the model involved the M-theory duality on 
AdSyxS'', Witten argued that it could equivalently be described as the background of a stack 
of Nc D4 branes in type IIA supergravity, where one of the directions parallel to the D4 
branes is compactified into a circle. Antiperiodic boundary conditions around the circle are 
imposed on the fermionic fields, giving the fermions a mass and breaking the supersymmetry. 
The scalar fields also acquire a mass at the one-loop level, leaving the SU{Nc) vector field 
as the only massless field in the theory, and thus reproducing nonsupersymmetric pure 
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Yang- Mills theory at energies small compared to the Kaluza-Klein scale Mkk- 

Sakai and Sugimoto added massless flavor to the theory by considering the addition of 
Nf D8-D8 probe branes to the background 1^, where the probe branes are transverse to 
the circle of compactiflcation. The essential idea of adding flavor to holographic systems via 
probe branes ISj is that if Nf <^ Nc, the backreaction of the probe branes on the geometry 
can be neglected, and the probe brane action consists simply of the Dirac-Born-Infeld (DBI) 
action in the original background, plus the relevant Chern-Simons (CS) terms. 

There is a U(A^) xU(A'^) gauge symmetry living on the probe brane pairs, which pro- 
vides a global chiral symmetry on the fleld theory side. The geometry of the D4 branes is 
cigar-shaped, and one flnds from analysis of the DBI action that the D8-D8 branes merge 
at some value of the radial coordinate. Thus, the l]{Nf)xl][Nf) symmetry is broken to a 
single \J(Nf)\ this is a holographic manifestation of chiral symmetry breaking. It was shown 
in Refs. [l6| that the DBI action of the probe branes reproduces much of the low-energy 
physics of QCD, including the chiral Lagrangian and qualitative features of the meson spec- 
trum. Inclusion of the CS term for the probe branes incorporates the effects of the axial 
anomaly into the low energy theory. 



It is possible to construct baryons in this model [17|. Witten's baryon vertex appears as 
a D4 brane wrapping the S^. Strings stretching between the D4 brane and the D8 branes 
will source the gauge fleld living on the D8s. The baryon number is then given in terms of 
the SU{Nf) valued fleld strengths F living on the probe branes as 

= tJ-t / d^xdz eMNpqtT [FmnFpq] , (1) 
SZn^ J 

where M, N,P,Q = 1, 2, 3, z. The baryon number integral also shows up in the probe brane 
action in the CS term coupled to the time-component of the U{1) part of the gauge fleld. A: 

Scs ^ -^^(^MNPQ / d'^xdz Aotr [Fmn^pq] . (2) 

Thus, Aq acts as a source for baryon number, and, in order to turn on a flnite chemical 
potential, we will consider solutions with nontrivial Aq. 

We also wish to turn on an external magnetic fleld that couples to our flavor degrees of 
freedom. There is no proper U(l) gauge fleld in our theory. However, we can simulate the 
effects of an external fleld by weakly gauging U(l) subgroups of the global chiral symmetry. 
In real QCD with Nf = 2 flavors, the electric charge is related to both the third component of 



isospin and to the baryon number, Q = + \B. Therefore, in order to properly introduce 
an electromagnetic field into Sakai-Sugimoto, one would have to gauge both one of the 
components of /a and the baryon number potential. However, for simplicity, in this paper 
we will only consider gauging of the isospin component by looking for solutions with a 
nonzero F^'^\ 



III. SOLUTION WITH FINITE MAGNETIC FIELD AND CHEMICAL POTEN- 
TIAL 

We will first establish our conventions and notation. 

In this paper, we will consider the case of two fiavor D8 branes. The branes and the 
antibranes will be maximally separated around the circle of compactification, as in the 
original treatment by Sakai and Sugimoto j^. 

The effective action of the probe branes is written in terms of a U(2) = U(1)bX SU(2) 
five-dimensional gauge field. It will often be convenient to distinguish the fifth spacetime 
index z from the four- dimensional boundary indices. Thus, upper-case roman letters M, A^... 
run over all five spacetime directions, whereas lower-case greek letters fi,^... run only over 
0, 1, 2, 3. We will use a combination of form and component notation to describe the gauge 
fields. The U(2)-valued form fields are given by 

A = AAidx^^, (3) 
J^ = dA + iAAA. (4) 

The gauge fields can be decomposed into the U{1) part A and the SU{2) part A as 

A = A + ^ = A^^r + ^1, (5) 

where the T* are the SU (2) generators normalized as tr \T^T^] = ^6''^ and the numerical 
factor on the U{1) piece is to ensure that the U{1) generator is normalized in the same 
manner. Likewise, is decomposed in terms of F and F: 

T = F^^T + ^. (6) 
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In component form, 



_ 4 (*) _ f),, A (») _ Jjk A (i) A 



MN 



SmAm - SmA 



The effective theory of the branes is described by the action ^ 



S — Sym + Scs, 
Sym = —K, I d^xdzii 



s 



cs 



2477^ 



tr 



where 



K 



2167r3 



and h{z) and k{z) are defined as 



hiz) = {1 + z^)-^^^ k{z) = l + z''. 



(7) 
(8) 

(9a) 
(9b) 

(9c) 
(10) 

(11) 



In "bottom-up" AdS/QCD models where the chiral symmetry is spontaneously broken by the 
boundary conditions at the IR brane, the action has the same form as (El), but with different 
functions h(z) and k{z). For example, in the model considered in [10|, h{z) = const and 
g{z) ~ cosh (62;), where 6 is a constant. 

The Yang-Mills action (19bp arises from expanding the DBI action for the probe branes 
to second order in field strengths. For large values of the field strength, this expansion will 
no be longer valid. We will see later that turning on a magnetic involves setting f1^2^ = B. 
It can be shown that the cubic terms involving B can be dropped as long as B satisfies the 
inequality 



27n B 



< 1. 



(12) 



Alternatively, the action (19bp can be interpreted as a bottom-up effective action. From this 
perspective, B is allowed to be arbitrarily large. 



In terms of the U{1) and SU{2) pieces, the action reads [l7 | 

h{z 



Sym 
Scs 



-— I d'^xdz 
2 



247r2 



3 - 1 1 

-A tiF^ + -AF"^ + -d{ Ail 
2 4 2 



2FA - -A^ 



(13) 
(14) 



In 171] , a localized soliton solution is found whose size is C>(l/v^). In this case, the scaling; 
of the gauge field solution allows the equations of motion to be expanded as a series in A. 
In our case, there is no such scaling in our solution, so we must work with the full equations 
of motion. Variation of the above action leads to the equations 

-Ml^k{z)dMF''-^ = -e^^f^\Fl^)F^l + (16) 
a 

^ [h{z)DMF^^^''' + Mi^9.(M^)F«^^)] = -e^^'-^^F.pFW, (17) 

^Mlj,k{z)DMF^^^^ = -e^'^fx^^F^pF^l (18) 

where the covariant derivative acting on field strengths is 

U^f^ii)^-^ = df.F'^^^" + e'^'' F^'''^^" . (19) 

and we have defined, for future convenience, 

Nc 277r 

lOTT^K 2A 

We wish to turn on a magnetic field and a baryon number chemical potential and to look 
for a solution homogeneous in Minkowski space. As discussed above, in order to turn on the 
magentic field we will assume a nonzero , generate a baryon number chemical 

potential we will assume a nonzero Aq. Recall that in terms of the five-dimensional gauge 
fields, the baryon number is given by ([1]), 

iV^ = ^ / d'x dz 6^'™Fi?^i^^ (21) 

where M, A^, P, Q are not zero. We see that in order to have nonzero baryon number, if FI2 
is nonzero then F^J must also be nonzero. Because we want a solution that is homogeneous 
in the four-dimensional coordinates, all these quantities will only depend on z. Let us assume 
all other fields are zero. To reduce clutter, we will henceforth denote F'-^-* simply by F. 

The second and fourth equations of motion are then trivially satisfied, whereas the first 
and third become 

^hLd,{k{z)F'^) = -F'^F^' (22) 
^{d,{k{z)F^^) = (23) 
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The covariant derivative has reduced to a partial derivative because is zero for j = 1, 2. 
We notice that due to the Bianchi identity and the requirement that fields depend only on 
Fi2 has to be a constant, F12 = B. The equations can be solved exactly for any function 

k{z). For k{z) = 1 + 2;^ as in the Sakai-Sugimoto model, the general solution is 



1 + 1 + 

exp(i? arctan^;) ^ exp(— arctan 2;) ^^^-^ 



where B = aB/Mj^^^. 

In the presence of a finite chemical potential, the thermodynamic ground state will min- 
imize the free energy H — ^Nb- Minimizing this quantity will give us the values of ci and 
C2 in the ground state. Under our Ansatz, the baryon number ([1]) reduces to 

= A / = T^^ix(ci + C2) sinh ^ (26) 

where V is the volume of the three dimensional space. 
The energy of the configuration is given by 



^ j d'xdz (^h{z)B' + Ml^k{z) (f3,F3^ - F,oF^o)^ . (27) 



The piece proportional to B gives an infinite contribution. This is the expected divergent 
energy of a space-filling magnetic field. This piece is independent of the constants ci and 
C2, so does not affect our minimization problem. 

Performing the integrals involved in the energy we write 

9 , 7iB I 2k , TiB ,99, u , ,1 , , 

H-^^NB = VMIj, sinh — I cosh —{cj + 4) - ^(ci + c^) j . (28) 

It is simple to minimize this with respect to ci and C2. The solution is 

ci = C2 = — (29) 

167r2/tM|^ cosh ^ 

Plugging into fl24p and fl25l) . we can write the ground state solutions entirely in terms of the 
parameters of the problem, 



p3z 



277r fiB ^Q^h(^^arctanz) 
^^i^ cosh (^4-) (1 + .2) 



(30) 



and 

277r fiB sinh(^^arctan2;) 



(31) 



The energy can also be found in terms of fundamental quantities by plugging in the values 
of Ci and C2 into our earlier expression for H. Doing so gives 

e = T7 = ^tanh --^ . (32) 



V 4:71^ Nc V 4A M|.^y 
One can express the energy density in terms of the baryon number density ub, 

nl , /277r2 B 

72 

The asymptotics of (133|1 at small B is 



e^.^iV.^coth(^^). (33) 



^ 4AiV, nlMl^ 

27 52 ■ ^ 



and at large B is 



2 



e 



vr^iV,:^. (35) 

Note, however, that in order to obtain the large B asymptotics we must assume that 
B/ {XM]^j^) ^ 1. This is precisely the limit in which we can no longer trust the quadratic 
approximation to the DBI action. Thus, in order to consider this limit, we must be think- 
ing in the context of a bottom-up AdS/QCD model, where higher order terms in F are 
suppressed from the beginning. 

IV. INTERPRETATION OF RESULTS 

We now try to interpret our results for small and large magnetic field B. At small 
magnetic fields, one can use chiral perturbation theory. The Hamiltonian describing the 



interaction of the vr*^ field with the magnetic field and baryon chemical potential is \l] 



(36) 



The minimum of H' is achieved at 



V.»^^B (37) 
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at which the energy and baryon number densities are 



The relationship between e and is 



nB = ,ro - (38) 



e = 8-V.^§. (39) 
Now using the value for the pion decay constant found in {4], 

= ^.Kk, (40) 

we reproduce the low-i? asymptotics of Eq. (!33l) exactly. 

We now show that at very large B the thermodynamic relation between e and ns, obtained 
in the approximation where one replaces the DBI action by the Maxwell action, approaches 
that of a free noninteracting gas. Consider a system of free, noninteracting quarks of two 
flavors u and d in an external magnetic field coupled to isospin. The charges of the quarks 
are 1/2 and —1/2. In a magnetic field, the transverse motion is quantized. The lowest 
Landau level consist of two branches, 

E = (41) 

each having degeneracy in the transverse direction equal to \e\B/{2'K) = B/{A7i). 

A baryon chemical potential ^ corresponds to a quark chemical potential fig = fi/Nc- The 
quarks fill energy levels below fig, leading to a nonzero baryon densty. At small chemical 
potentials (or, equivalently, large magnetic field) the filled energy levels are all in the lowest 
Landau level. The total baryon number is then 

B f dk, fiB 



-fJ,q 



Inverting the relation. 



/i = ^ ub (43) 



and integrating over hb, one finds the energy density as a function of the baryon number 
density 



n'N^^. (44) 
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This coincides exactly with the thermodynamics of our bottom-up model at large B. There- 
fore, we conclude that at large B, the equation of state of matter at finite-baryon density is 
identical to that of free quarks. 

Moreover, by redoing the previous calculation, one can show that this feature is indepen- 
dent of the choice of the function k(z), given that the integral f^dkk~^{z) is convergent. 



This condition is satified in the model of 



10|] where k{z) ~ cosh.{bz). Therefore, in bottom-up 



holographic models of QCD where the action contains only the Maxwell and Chern-Simons 
parts, and where chiral symmetry breaking is due to a boundary condition at 2; = 0, the 
equation of state at very high magnetic field approaches that of a free gas. 



V. CONCLUSIONS 



To conclude, we have found a solution to the field equations in the Sakai-Sugimoto model 
that corresponds to matter at finite baryon density in an external magnetic field. In contrast 
to the case without a magnetic field, it is possible to write down a solution that is completely 
homogeneous in space. At small B and small chemical potential the result can be obtained 
from the chiral Lagrangian with the Chern-Simons term. The solution continues to exist 
for any value of B, although for larger values of B the solution can no longer be viewed as 
arising from the full Sakai-Sugimoto model. We can, however, interpret our solution at large 
B as arising from a bottom-up AdS/QCD model. At large B the system behaves, from the 
point of view of zero-temperature thermodynamics, as a system of free quarks. 

How one can explain the latter fact? Right now we have only some vague ideas of how 
it can be understood. In magnetic fields the fermions move in Larmor orbits whose radius 
shrinks as B —>■ 00. This fact may explain why interaction between quarks do not seem to 
play any role at large B. 

Clearly, the situation should be investigated further. One question one can ask is whether 
the state is stable with respect to small perturbations. At small /i, the energy per baryon is 
small and the system is clearly more stable than the ordinary Skyrmion crystal. We leave 
the investigation of the stability of the system at large B and fi to future work. 
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